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ABSTRACT 

Timing noise in the data on accretion-powered millisecond pulsars (AMP) appears as irregular pulse phase 
jumps on timescales from hours to weeks. A large systematic phase drift is also observed in the first discovered 
AMP SAX J1808.4— 3658. To study the origin of these timing features, we use here the data of the well 
studied 2002 outburst of SAX J1808.4— 3658. We develop first a model for pulse profile formation accounting 
for the screening of the antipodal emitting spot by the accretion disk. We demonstrate that the variations of the 
visibility of the antipodal spot associated with the receding accretion disk cause a systematic shift in Fourier 
phases, observed together with the changes in the pulse form. We show that a strong secondary maximum 
can be observed only in a narrow intervals of inner disk radii, which explains the very short appearance of 
the double-peaked profiles in SAX 11808.4—3658. By directly fitting the pulse profile shapes with our model, 
we find that the main parameters of the emitting spot such as its mean latitude and longitude as well as the 
emissivity pattern change irregularly causing small shifts in pulse phase, and the strong profile variations are 
caused by the increasing inner disk radius. We finally notice that significant variations in the pulse profiles in 
the 2002 and 2008 outbursts of SAX J1808.4— 3658 happen at fluxes differing by a factor of 2, which can be 
explained if the inner disk radius is not a simple function of the accretion rate, but depends on the previous 
history. 

Subject headings: accretion, accretion disks - methods: data analysis - pulsars: individual (SAX 
J1808.4— 3658) - stars: neutron - X-rays: binaries 



1. INTRODUCTION 

The first accreting millisecond X-ray pulsar (here- 
after AMP) SAX J1808 4 -3658 w as discovered in 1998 
dWiinands & van der Khsl 119981: IChakrabart v & Morpaiil 
119981) . Since then, 12 objects of this class have been dis- 
covered, with the spin frequencies in the range 182-599 Hz. 
Many AMPs (e.g., IGR J00291H-5934, XTE J1751-205, 
XTE J1807-294, SAX J1808.4-3658, XTE J1814-338) 
often show nearly sinusoidal profiles, consistent with only 
one visible hotspot, probably because the accretion disk 
extends very close to the stellar surface. The general stability 
of the pulse profile allows to obtain a high photon statistics 
and to use the average pulse profile to get constraints on the 
neutron star mass-radius relation and the equation of state 
dPoutanen & Gierliriski 2003). 

However, sometimes the profiles do show variations. In 
SAX J1808.4— 3658, for example, a dip appears around the 
pulse maximum at high photon fluxes, and at low fluxes the 
profile becomes skewed to the left and sometimes is double- 
peaked. The pulse profile evolution during its five observed 
outbursts is amazingly similar, with significant changes in the 
profiles accompanied by large (0.2 cycles) jumps in the phase 
of the fundamental (see Burderi et al. 2006; Hartman et al. 
20081 120091; iPatruno et al.i i2009ai: ilbragimov & Poutaneni 



20091 hereafter IP09). 



Besides these noticeable profile changes, the timing noise 
with timescales from hours to weeks is known to be present 
in the data of this and other AMPs, when the measured 
pulse phase delays s how fluctuations around a mean trend 
jPapitto et al. 2007; IChou et al.1 llOOSt iRiggio et all l200l 
IPatruno et al.n2009h) . There are many possible reasons for 
the pulse profile changes and the timing noise: changes in the 



spot shape and position, variations of the emission pattern, in- 
ner disk radius o r the optical depth of the accretion stream 
( Poutaneni 120081). Changes in the posi tion of the hotspots 
jLamb et al. 1 120081 IPatruno et al.ll2009bl) can cause random 
phase irregularities; however, it is difficult to imagine that 
they can lead to qualitative changes in the pulse form (such as 
transformation of a single-peaked profile to a double-peaked) 
unless accompanied by strong variation in the emissivity pat- 
tern. 

On the other hand, the visibility of the secondary spot might 
change even with small fluctuations of the accretion rate. IP09 
proposed that the strong pulse profile variations in the 2002 
outburst of SAX J1808.4— 3658 results from the appearance 
of the secondary antipodal spot to our view when the magne- 
tospheric radius increases and the accretion disk recedes from 
the neutron star with the dropping accretion rate. The decreas- 
ing amplitude of Compton reflection in the spectra during the 
outburst supports this interpretation. 

In this Letter, we directly fit the pulse profiles of SAX 
J1808.4-3658 during its 2002 outburst with the theoretical 
model for AMP pulses including the effect of partial screen- 
ing of the antipodal emitting spot by the accretion disk. We 
aim at quantifying the variations in the position of the hotspot 
centroid, its emissivity pattern and the inner accretion disk ra- 
dius, and determining the causes of the variations in the pulse 
shape and the timing noise. 

2. MODEL 

We assume that the stellar magnetic dipole is displaced by 
angle 6 from the rotational axis. The accretion disk is dis- 
rupted by the stellar magnetic field at some radius which 
is a function of the accretion rate and the neutron star mag- 
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Figure 1. Pulse profiles produced by two antipodal spots of vaiious shapes 
for iJin varying from 20 to 30 km with step 2.5 km (from bottom to top). The 
left hand side panels con'espond to the blackbody emissivity pattern h = 0, 
and the panels at the right are for h = 0.8. At lowest radii, the view of 
antipodal spot is almost fully blocked by the disk, while at highest radii, the 
whole spot is visible. Stellar parameters: M = IAMq, R = 10.3 km, 
i = 65°, 9 = 10°. The spot angular size p is given by equation (T}, and for 
the ring and crescent geometries, the inner radius is p/ y/2. 

netic moment. Matter following magnetic field lines hits the 
stellar surface making two antipodal hot spots. For the dipole 
field, the outer spot edge is displaced from the magnetic pole 
by an angle p: 



sm p = cos 

h 



(1) 



where R is the neutron star radius. This gives us a rough esti- 
mation of the spot size as a function of i^in. The spot size and 
shape depend on the magnetic inclination and the inner disk 
radius in a complicated fashion. When 0^1 and » R, 
the hotspot is nearly a circle. When the disk is close to the 
star, the spot becomes a ring-like because the matter accretes 
via field lines penetrating the disk within a narrow ring. If 
the magnetic di pole is sufficie ntly inclin ed, the spot becomes 
a cres cent-like ( Romanova et al . 2004; Kulkarni & Romanoval 
120051) . Thus we consider three simple geometries which in- 
corporate all these cases: (1) a homogeneously emitting cir- 
cle, (2) a ring, and (3) a crescent, which we approximate as a 
part of the ring. 

For a given spot geometry and stellar parameters, we then 
compute the pulse profile. We account for general and special 
relativity effects (such as gravitational light bending, redshift, 
aberration, Doppler boosting), and time delays. We also need 
to specify the spectral shape and the angular emissivity pattern 
of radiation from the spot. We follow the general methodol- 
ogy as described in detail in Poutanen & Gierliiiski (2003), 
iViironen & PoutanenI ( 12004ft . and Poutanen & Beloborodoj 
(120061) . As the accretion disk might be rather close to the 
stellar surface, we account also for the absorption by the disk 
as described in Appendix C of IP09. We also approximate the 
theoretical profile by a sum of two harmonics: 

= F{l+aiCOs[27r(0-0i)]+a2COs[47r(0-02)]}- (2) 
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Figure 2. Fourier amplitudes and phases of the fundamental and overtone (as 
given by Equation ^2)) of the pulse profiles shown in Figure [T] as functions 
of the inner disk radius. Solid and dashed curves correspond to the cases of 
h = and h = 0.8, respectively. 

For the illustration, we assume, rather arbitrarily, the fol- 
lowing parameters: neutron star mass AI — IAMq, radius 
R = 10.3 km, spin frequency 401 Hz (corresponding to SAX 
J1808.4— 3658), inclination i = 65°, magnetic inclination 
9 ~ 10°, and a spot of angular radius p given by Equation ([T]). 
The spectral energy and angular distributions of the radiation 
intensity at the stellar surface are assumed to be represented 
by the expression 

I{E,a)(x {1-h cos a) E-^^-^\ (3) 

where F w 2 is the power-law photon spectral index observed 
in AMPs and SAX J1808.4-3658 in particular (see, e.g. 
Pouta nen & Gierliiiski 2003; Falangaetal. 2005a b, 200"| 
iPoutane n 2006), a is the angle relative to the surface normal, 
and h is the anisotropy parameter The assumed angular de- 
pendence with h = corresponds to the blackbody-like emis- 
sion pattern, while h « 0.6-0.8 is consistent with the thermal 
Comptonization in an optically thin slab, which is probably 
responsible for the power-law-like spectra. 

For a ring- and a crescent-like spots we assume the inner 
radius of the spot p/V^, and for the crescent we consider 
half a ring furthest away from the ro tational pole similar to 
what is obtained in MHD s imulations dRomanova et al.ll2004l ; 
iKulkarni & Romano vall2005i) . Figure [T] shows the pulse pro- 
file variations as a function of the inner disc radius for the 
three geometries and two anisotropy parameters h — 0, 0.8. 
The noticeable signatures of the second spot appear when 
i?in > 20 km, the dramatic change in the pulse profile hap- 
pens when the inner radius changes from 24 to 26 km, as 
most of the antipodal spot then appears to our view. These 
fast changes are also reflected in a rather large phase shift of 
the fundamental by ~0.2 (see Figure|2l), but generally smaller 
variations of the phase of the overtone, similar to what is ob- 
served in SAX J1808.4-3658 (see Figures 6 and 7 in IP09). 
The behavior of the amplitudes depends on the anisotropy. 

While blackbody pulse shapes are rather smooth, the 
Comptonization emissivity pattern lead to the presence of nar- 
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Figures. Sample of pulse profiles of SAX J18()8.4— 3658 observed during the October-November 2002 outburst, in order of decreasing flux. The top and 
bottom panels coiTespond to the profiles in the 3.7-5.7 keV and 9.8-23.2 keV energy bands, respectively. The coiTesponding best-fit models parameters are 
indicated. The solid curves give the modeled pulse shapes. 



rower features and more pronounced secondary maxima; the 
ring generates sharper features in the lightcurves and more 
rapid evolution with i?in comparing to the circle case. The 
crescent gives higher variability amplitude, because of its 
larger mean displacement with respect to the rotational pole 
of the star. However, the initial (when only one spot is visi- 
ble) and the final (when both spots are visible) pulse profiles 
depend very little on the spot geometry, and changes of the 
Fourier phases with inner radius are also very similar We 
also note that clearly double-peaked profiles are produced in 
a rather narrow range of when only a fraction of the sec- 
ondary spot is visible. At large i?i„, the profiles becomes 
again nearly sinusoidal in spite of the fact we observe both 
spots simultaneously. 

The phase shift of the fundamental shown in Figure |2] is a 
very robust prediction independent of the stellar parameters, 
emissivity pattern, and spot shape. The exact disk radius when 
this occurs depends on the stellar compactness and observer 
inclination (see 1P09). For a sufficiently small i and a not very 
compact star, the antipodal spot might never be visible. In that 
case, there will be no noticeable changes in the pulse profile 
with the accretion rate. 

3. COMPARISON WITH THE OBSERVED PROFILE EVOLUTION 

Now we make a detailed comparison of the theoretical pulse 
profiles predicted by the model to the data on 2002 outburst of 
SAX J1808.4-3658 well covered by the Rossi X-ray Timing 
Explorer. The data and their reduction are described in details 
in 1P09. At the peak stage (see Haitman et al. 2008 and 1P09 
for the definition of outburst stages), profiles are nearly sinu- 
soidal with a dip noticeable at the pulse maximum at energies 
above 10 keV. This might be caused by absorption in the ac- 
cretion column (see IP09) complicating the pulse modeling, 
and therefore we do not consider this stage in the paper. Dur- 
ing the slow decay stage on October 17-27, the pulse profile, 
slightly skewed to the right (see Figure[3] left panels), was ex- 
tremely stable, in spite of the fact that the flux was changing 
by a factor of 4. At the rapid drop stage (October 27 - Novem- 
ber 1) the flux fell down rapidly, while the pulse profile was 
evolving: on October 29 the secondary maximum - a signa- 
ture of the second hot spot - becomes pronounced and after 



October 30 the pulse became skewed to the left and the sec- 
ondary maximum shifted to the rising part of the main peak. 
In Figure[3] we present a subset of the observed pulse shapes. 

The pulse profiles can be rather well described by a sum of 
two Fourier harmonics, fundamental and overtone, given by 
Equation Q. The dramatic variability of the pulse profiles is 
reflected in a phase shift of the fundamental by 0.2 observed 
within a couple of days around October 30, whil e the phase 
of the overtone did not change significantly ( Burderi et alj 
2006; Hartman e t al. 2008; 1P09). The behavior of the pulse 
phases and amplitudes varies with energy: for example, 
the amplitude of the high-energy pulse decreases with time, 
while the pulsation amplitude at low energies noticeably in- 
creases (1P09). The reason for that most probably lies in the 
fact that the soft blackbody which c ontributes significantly 
to the total spectrum below ^7 keV dGierliriski et al.l 120021 ; 
lPoutaneir& Gierlinski 2003; 1P09) might be produced in a re- 
gion displaced from the accretion shock where harder X-rays 
are produced. 

For comparison of the model to the data, we take the same 
stellar parameters as given in Section |2l inclination i = 60°, 
and consider a filled circle geometry of the spot. We fit simul- 
taneously the pulse profiles in two broad energy bands 3.7- 
5.7 and 9.8-23.2 keV. The second band is dominated by the 
Comptonized components, while the blackbody contributes at 
the lower-energy band. We consider a possibility that photons 
at these energies can be produced predominantly at different 
positions, determined by the colatitude Og.h and the azimuth 
0s, h (in units of 27r) of the spot center The angular emissiv- 
ities of the radiation in the two bands are described by Equa- 
tion ^ with parameters /is,h- The inner disk radius i?in is of 
course the same for both energies at each moment of the out- 
burst, and the spot size is given by Equation ([T]) with = 0. 
The best-fit parameters are given in Table [T] and the represen- 
tative pulse profiles are shown in Figure [3] 

We see that the model describes the pulse profiles rather 
well, except some narrow features. In general, the resulting 
i?in increases when the flux decreases, but not monotonically. 
We also see some variations in positions of the spot centroids, 
both in latitude and in longitude. These changes can be as- 
sociated with the timing noise (i.e., erratic variations of the 
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Table 1 

Best Fit Parameters 



MJD 


Date 


F" 








<t>h 


hs 




Rin 


xVd.o.f. 








(deg) 


(deg) 










(km) 




52569.0-52569.5 


October 22 


15.6 


10.2 


10.8 


0.10 


0.13 


0.68 


0.80 


16.8 


34/56 


52580.8-52581.5 


November 3 


4.2 


9.7 


5.2 


0.14 


0.10 


0.79 


0.85 


21.7 


80/56 


52581.6-52582.0 


November 3 


3.6 


8.2 


5.4 


0.15 


0.06 


0.75 


0.83 


22.5 


62/56 


52576.9-52577.4 


October 30 


3.5 


11.6 


8.1 


0.13 


0.08 


0.72 


0.78 


21.3 


86/56 


52577.9-52578.1 


October 3 1 


2.4 


18.4 


6.6 


0.13 


0.09 


0.67 


0.84 


22.8 


92/56 


52583.8-52584.2 


November 6 


1.9 


18.1 


7.0 


0.15 


0.12 


0.68 


0.82 


22.7 


61/56 


52578.9-52579.5 


November 1 


1.5 


19.3 


5.7 


0.12 


0.11 


0.61 


0.79 


23.2 


74/56 



' Flux in 2.5-25 keV band in units 10 



erg cm 



Fourier phases). The jump in the fundamental phase observed 
around October 30, which happened together with the dra- 
matic change in the pulse profile, is explained here by the 
appearance of the antipodal spot when the inner disk radius 
is about 20 km (the exact value would depend on the stellar 
compactness and the inclination). The origin of some dis- 
crepancies between the model and the data probably lies in a 
simplified spot shape assumed here (the pulse profiles depend 
considerably on that, see Section |2]i and in the assumption of 
the equal area spots producing emission at different energies. 

Since the evolution of the pulse profiles during various out- 
bursts of SAX J1808.4— 3658 is very similar, our model thus 
can be applied directly to other outbursts too. It is worth 
noting, however, that the changes in the pulse profile and 
appearance of double-peaked profile in the 2008 data hap- 
pen at a flux twice as large as in the 2002 and 2005 out- 
bursts (see iHartman et al.l 120081 [2009). This might imply 
that the inner disk radius and the magnetospheric radius are 
not a simple function of the accretion rate, but also of the 
history of the outburst. We can speculate that in 2008, a 
smaller total accreted mass resulted in a less effective re- 
duction of the stellar magnetic field (Gumming 2008) and, 
therefore, the same inner disk radius could have been reached 
at a higher accretion rate. The phase fluctuations correlated 
with the flux as obser ved, for exarnple, in XTE J 1807— 294 
andXTEJl 8 14-338 ('Papitto et al. 20071: laiouetaLl 120081: 
iRiggio etal. 2008: Patruno et al. 2009bJ, cannot be easily ex- 
plained by the inner disk variations. 

4. CONCLUSIONS 

We have developed a model for the pulse profile formation 
in AMPs which accounts for the partial screening of the an- 
tipodal emitting spot by the accretion disk. We have demon- 
strated that the appearance of the antipodal spot, due to the 
increasing inner disk radius, leads to sharp changes in the 
pulse profile and corresponding jump in Fourier phases. We 
showed that a strong secondary maximum appears only in 
a rather narrow intervals of inner disk radii, explaining the 
very short appearance of the double-peaked profiles in SAX 
J1808.4-3658. 

By directly fitting the model to the pulse shapes of SAX 
J1808.4— 3658 observed during its 2002 outburst, we were 
able to quantify the variations in the position of the spot cen- 
troids, their emissivity pattern and the inner accretion disk ra- 
dius. The sharp jumps in Fourier phases, observed together 
with the strong changes in the pulse form, can be explained 
by the variations of the visibility of the antipodal spot asso- 
ciated with the receding accretion disk. We also note that a 
factor of 2 difference in fluxes in 2002 and 2008 outbursts of 



SAX J1808.4-3658, when the pulse profile started to change 
significantly, implies that the inner disk radius is not a simple 
function of the accretion rate, but might depend on the previ- 
ous history. We finally note that many physical timing noise 
mechanisms probably operate in AMPs simultaneously. 
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